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ABSTRACT 

We present a detailed study of the QSO-galaxy pair [SDSS J163956. 35+112758.7 (z q = 
0.993) and SDSS J163956.38+112802.1 (z g = 0.079)] based on observations carried 
out using the Giant Meterwave Radio Telescope (GMRT), the Very Large Baseline 
Array (VLB A), the Sloan Digital Sky Survey (SDSS) and the ESO New Technology 
Telescope (NTT). We show that the interstellar medium of the galaxy probed by the 
QSO line of sight has near-solar metallicity (12+log(0/H) = 8.47±0.25) and dust 
extinction (E(B-V)~ 0.83±0.11) typical of what is usually seen in translucent clouds. 
We report the detection of absorption in the A6284 diffuse interstellar band (DIB) 
with a rest equivalent width of 1.45±0.20A. Our GMRT spectrum shows a strong 
21-cm absorption at the redshift of the galaxy with an integrated optical depth of 
15.70±0.13 km s _1 . Follow-up VLBA observations show that the background radio 
source is resolved into three components with a maximum projected separation of 
89 pc at the redshift of the galaxy. One of these components is too weak to provide 
useful 21-cm H I information. The integrated H I optical depth towards the other 
two components are higher than that measured in our GMRT spectrum and differ 
by a factor 2. By comparing the GMRT and VLBA spectra we show the presence of 
structures in the 21-cm optical depth on parsec scales. We discuss the implications of 
such structures for the spin-temperature measurements in high-z damped Lyman-a 
systems. The analysis presented here suggests that this QSO-galaxy pair is an ideal 
target for studying the DIBs and molecular species using future observations in optical 
and radio wavebands. 

Key words: quasars: active - quasars: absorption lines - quasars: individual: SDSS 
J163956.35+112758.7 - ISM: lines and bands - ISM: molecules 



1 INTRODUCTION 

Understanding parsec scale H I opacity fluctuations in the 
interstellar medium (ISM) of galaxies and how they de- 
pend on the feedback from in situ star formation is very 
important. In the Galaxy, the radio observations of H I 21- 
cm absorption towards high-velocity pulsars and extended 
radio sources, the optical observations of Na I absorption 
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lines towards globular clusters and binary stars have shown 
that the diffuse ISM is structured on parsec to sub - parsec 
(~10AU) scales (IFrail et al.lll994l; IBates et aT]|l995l: iHeileU 
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Deshpandd |2000|; iRollinde et al.1 [ 2003: 



Brogan et al 



1200a 



Roy et al 



20121 ). 



Boisse et alj 



These small- 
scale structures, due to large over-pressures, cannot sur- 
vive rapid evaporation in the s tandard pressure-equilib rium 
based models (see for example. iMcKee fc Ostrikerlll977l ) and 
therefore raises many important questions regarding the na- 
ture of the ISM. It is possi ble that these are transient phases 
controlled by turbulence (|Mac Low fc K lcsscn 20041) . 

Physical conditions in H I gas associated with high red- 
shift galaxies can be probed by 21-cm absorption. In par- 
ticular, the H I 21-cm optical depth combined with a deter- 
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mination of iV(H i) from Lyman— a absorption yield an es- 
timat e of the spin temperature, Ts, (|Kanekar fc Chengalurl 
l2003h which is a good indicator of the kinetic temperature 
l|Rov et all 1200(1 . However, to fully interpret high redshift 
21-cm spin-temperature measurements one needs to know 
the parsec scale structure of the H I gas. In turn, if radio 
emission of the background QSO has structures at differ- 
ent scales then Very Large Baseline Interferometry (VLBI) 
spectroscopy can be used to probe the spatial variations of 
the H I gas opacity. 

Due to limitations of existing facilities, VLBI spectro- 
scopic observations of QSO absorbers are limited to low- 
redshifts onl y. In the ca s e of t he z a bs= 0.0912 DLA towards 
B0738+313, lLane et all (|2000h found the background source 
to be partially resolved at milliarcsec (mas) scales. Within 
the measurement uncertainties they do not find any strong 
variations in the H I optical depth across 20 pc. The asso- 
ciated galaxy could be a low surface brightness dwarf at an 
impact parameter of ^3. 5 kpc without any sig nature of on- 
going star formation (see iTurnshek et al.ll2001 ). In the case 
of the 2 = 0.03321 galaxy towards J104257.58+074850.5, 
the background radio source is unresolved in the Very Large 
Baseline Array (VLB A) image and the 21- cm absorption is 
found to be similar over 27.1x13.9 pc (see iBorthakur et al.l 
|2010| . for details). Note that the impact parameter in this 
case is 1.7 kpc and the foreground galaxy is a low luminos- 
ity spiral with low gas phase metallicity and dust extinction 
along the QSO sight line. 

Alternatively, time variability of 21-cm absorption 
is seen in two ca ses [the absorber s towards B0235+164 
(I Wolfe et all 1 19821 ) and B1127-145 (|Kanekar fc Chengalurl 
l200ll )] suggesting that the H I gas is patchy on pc scales. 
Interestingly, the z a b s = 0.524 absorption system towards 
B0235+164 is one of the extra-ga lactic systems that show s 
Diffuse Interstellar Bands (DIBs; [junkkarinen et al]l2004l ). 
Using a large sample of z ^ 1.5 Mg n absorbers searched 
for 21-cm ab s orptio n together with mas scale radio images 
iGupta et al.1 |2012l ). have concluded that the H I gas is 
patchy with a typical correlation length of 30—100 pc based 
on different correlation analyses. 

It has been shown that, like 21-cm absorption, one 
can use H2 absorption to probe the cold neutral medium 
(CNM) in z ^ 1 . 8 damped Lyman— a systems (DLAs; 
Ledoux et al.ll2003l;ISrianand et"aL Il2005l ; I Noterdaeme et al.l 



Table 1. Emission line 
SDSS J163956.38+112802.1 



parameters of the galaxy 



2008). Recently iBalashev et al.l (|201l|), using partial cover 



age arguments, derived the linear size of the liVbearing core 
and the H I envelope to be 0.15^0 05 pc and 8.2I4 J pc re- 
spectively for the z ab B = 2.3377 DLA towards Q1232+082. 
ISrianand et~aH j2012tl , combining 21-cm and H2 information 
for a sample of high- 2 DLAs, concluded that the typical size 
of H2 bearing clouds is ^ 15 pc. 

All this indicates that the ISM of high-z galaxies is 
structured on parsec scales. However, most of the evidences 
are based on indirect arguments. In addition one would like 
to connect the structures seen in H I to other observable 
properties of the galaxy such as metallicity, dust content, 
local star formation etc. We are performing a systematic 
search for 21-cm and OH absorption in a well selected sam- 
ple of Quasar-galaxy pairs (QGP, a fortuitous alignment of a 
foreground visible galaxy with a distant background quasar) 
with impact parameter b ^ 20 kpc. In case of detection 
we also perform optical long-slit and VLBA spectroscopic 
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follow-up observations. Initial results from the GMRT ob- 
servations of 5 QGPs, with the redshift of th e galaxies in 
the ra nge 0.03^ z g ^ 0.18 are presented in I Gupta et al.l 
J2010l ). 

Here we present a detailed analysis of a very special 
QGP with the background QSO (J163956.35+112756.7, z q 
= 0.993) sight line piercing through the optical disk (Fig.[T} 
of the foreground galaxy (J163956. 38+112802.1, z 3 = 0.079 
or corresponding heliocentric velocity of Uhoiio ~ 23710 
km s _1 ). We report the detection of 21-cm, Na I absorp- 
tion and DIBs in this galaxy. We use a long slit spectrum 
to study the star formation rate, metallicity and rotation 
curve of the galaxy. The background radio source is well re- 
solved at mas scales. We present the VLBA spectroscopy 
and discuss the parsec scale structures in the H I gas. We 
will use a flat cosmolo gical model with fi m = 0.27, h = 0.70 
(|Komatsu et al.ll2009l L 



2 PHYSICAL CONDITIONS IN THE GALAXY 

2.1 Metallicity and dust in the emission line 
regions 

Strong emission lines from the foreground galaxy are seen in 
the SDSS spectrum (see Fig. [T]). They suggest the presence 
of a star forming region within a projected separation of 2.2 
kpc to the QSO sightlineQ. We fit the emission lines using 
Gaussians and the integrated line fluxes are summarised in 
Table [T] A luminosity distance of 354.1 Mpc is used to con- 
vert these fluxes into luminosities. The last column in this 
table gives the ratio of the Ha flux with respect to that of 
other emission lines. Following Argence & Lamareille (2009), 
we derive the optical depth in the systemic V-band of the 



galaxy, r v 



= 1.39 



+0.18 



using th e observed Hq/H/3 ra- 



tio, t he intrinsic Balmer ratio of 2.85 (|Osterbrock fc Ferlandl 
120061) and the wavelength depende nce of the dust op tical 
depth (i.e t\) as given by Eq. 3 of I Wild et ail (|2007m . We 
estimate a dust uncorrected surface density of the star for- 
mation rate (SFR) of 0.006 Meyr^kpc^ 2 . Applying the 
dust correction will increase this estimate by a factor of 



throe for 
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quoted above. If we assume that this 



The SDSS fiber has a diameter of 3 arcsec and was centered on 
the QSO. An angular scale of 1 arcsec corresponds to a projected 
size of 1.474 kpc at the redshift of the galaxy. Therefore, the SDSS 
spectrum is the juxtaposition of the QSO spectrum with that of 
the galaxy within a radius of 2.2 kpc. 
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Figure 1. Left: SDSS colour image (30" X 30") showing the QSO (SDSS J163956.35+112758.7, z q = 0.993) sight line piercing through 
the disk/spiral-arm of the foreground galaxy (SDSS J163956. 35+112802.1, z s =0.079) Right: The SDSS spectrum of the QSO is fitted 
with a combination of a reddened QSO composite spectrum and a spiral galaxy template (blue spectrum) from SDSS. The insert gives 
the distribution of Ay measured for a control sample of QSOs with similar redshift as Q1639+1127. The arrow marks the measured Ay 
in the present case. In the top right panel we show the relative contribution of the galaxy light to the SDSS spectrum as a function of 
wavelength. This contribution derived in different bands from SDSS images are also plotted as big black points. 



disk galaxy obeys the Kennicutt-Schmidt law then we get 
N(H i)^ 2 . 6 x 10 2 1 cm -2 using the best fitted parameters of 
iKennicuttl (jl998al fbh . Using the measured 03N2 ratio (i.e log 
[F A (0 iiiA5007)/F a (H/3) ~ F A ([N i i1A6585)/F a (Hq ) 1) and 
the best fitted relationship given in IPettini fc Pagell (|2004T ) 
we get, 12+log (O/H) = 8.47±0.25. Thus the QSO sightline 
is passing close to (i.e. within 2.2 kpc) a star forming region 
having near-solar metallicity and high iV(H i). 



2.2 Modelling the QSO spectral energy 
distribution 

The background QSO spectrum is highly reddened (see 
Fig. [T]). We get an independent estimate of ry and E(B- 
V) by fitting the QSO spectral energy distribution (SED) 
using, 



h = Af^e 



~ TA2 + Bfg 2 . 



Here fx , fg 1 and fg 2 are, respectively, the observed flux at 
A, the flux in the QSO composite spectrum at Ai = A(l + 2 g ) 
and the flux in the galaxy continuum that also enters the 
SDSS fiber. We approximate the latter with a spiral galaxy 
template at A2 = A(l + z g ). A and B are normalization 
factors and z q and z g are the redshifts of the QSO and 
the foreground galaxy respectively. T he form of t\ is taken 
from the Milky Way extinction curve (|Fitzpatrick fc Massal 
1 19881 ). The fitting method u sed is very similar to the on e 
used in lSrianand et all (|2008l ) and lNoterdaeme et all l|2009f ). 
As can be seen from the right panel of Fig. [TJ the best fit 
(with xi = 1-23) to the observed spectrum is obtained for 



A v = 2.56±0.02 (with E(B-V)= 0.82±0.01 and Ry = 3.10, 
where A\ — 1.086 t\. The quoted errors are mainly statis- 
tical ones. In Fig. [1] we also show the distribution of Ay 
obtained for a control sample of SDSS QSOs with z ~ z q . 
For this exercise we do not consider the foreground galaxy 
contribution to the SED (i.e the second term in Eq. 1). The 
rms of the Ay distribution of 0.11 reflects a typical sys- 
tematic error in the SED fitting method due to the disper- 
sion of the QSO unreddened spectral energy distribution. 
Therefore, the reddening noted in the present case is sig- 
nificant at more than the 20cr level. The measured value 
of Ay is consistent with the line of sight passing through 
a translucent region [defined as a region with 1^ A y ^10 
(Ivan Dishoeck fc Black! 119891 ; ISnow fc McCallll2006l )]. If we 
use the relatio nship between AT(H 1) and Ay derived in the 
Galactic ISM (|Bohlin et alj|l97gft then the above inferred 
Ay is consistent with N(R 1) = 4x 10 21 cm" 2 . The difference 
between the two Ay estimates obtained using the Balmer 
decrement (i.e A v almcr = 1.5lt^|) and the fit of the Q SO 
SED can be attributed to differences in the dust opacity 
(or differences in AT(H 1), if we assume uniform dust-to-gas 
ratio) towards the QSO and the Ha emitting regions. In 
passing, we wish to point out that the line of sight discussed 
here has more dust opacity than any of the high-z dusty 
Mg 11 and CO system s that also produce 2175 A absorption 
JSrianand et al. 120081 ; iNoterdaeme et aL 1 120091 ; IJiang et al.l 
12Q11J). 

The SED fitting procedure used above predicts the rel- 
ative contributions of the QSO and the galaxy as a function 



2 The derived value of E(B-V) depends on our choice of Ry. In 
the absence of Ry measurement we adopt Ry = 3.1 as seen in 
the Milky Way. 
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Figure 2. NTT spectrum of the QSO. Vertical dotted and dashed 
lines mark the positions of Na D lines and DIBs respectively. A 
strong absorption is seen at the expected position of the A6284 
DIB. The best gaussian profile fitted to this feature is shown to- 
gether with the abso rption in the A5780 A D IB expected from the 
correlation found bv lFriedman et alj j20 1 lT ) in the Galactic ISM. 
Top-panel: Residual flux after subtracting the Gaussian fits from 
the observed spectrum. The lcr range from the error spectrum 
(horizontal dashed line) is also shown. 



of wavelength (top right panel in Fig. Q]). We independently 
estimate this ratio at different wavelengths using different 
SDSS broad band images within an aperture of 1.5" cen- 
tered at the position of the quasar. We use the following 
steps: 1) cut an image stamp of 80 pixels on a side around 
the QGP, 2) subtract a constant background obtained by av- 
eraging the counts in more than 10 neighboring regions close 
to the pair, 3) obtain the radial surface brightness profile of 
the galaxy (after masking the quasar) by fitting isophotal 
ellipses using STSDAS package in IRAF, 4) derive the total 
counts due to the galaxy within an aperture of 1.5 arcsec ra- 
dius at the location of the quasar using the fitted isophots. 
The estimated ratios together with the errors (mainly from 
Poisson statistics) are also shown in the top-right panel of 
Fig. [1] They are very much consistent with the curve ob- 
tained from the SED fitting. This confirms the robustness 
of our fitting procedure. In addition, based on the ellipticity 
of the outer isophots in the r and i bands, we measure a 
disk inclination angle of ~ 63 deg. 

2.3 Na I absorption 

Using the European Southern Observatory (ESO) New 
Technology Telescope (NTT), we performed long-slit spec- 
troscopic observations of the QSO-galaxy pair aligning the 
slit (of width 1.2 arcsec and length 4.1 arcmin) along the 
semi-major axis of the galaxy. The spectra (4 x 2700s expo- 
sures) were obtained in good seeing conditions (~ 0.8 arc- 
sec) using the ESO Faint Object Spectrograph and Camera 
(EFOSC2) covering the wavelength range 6040-7135 A at 
a spectral resolution of 200—250 km s _1 . The data were 
processed using standard long-slit routines. One di- 

mensional spectra were extracted with sub-apertures having 

3 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which arc operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



a width of 5 pixels (i.e projected size of ~2 kpc). The results 
are summarized in Fig. [2] 

In the spectra of the QSO (shown in the bottom panel 
in Fig. [2]) we detect Na lAA5891,5897 lines with rest equiv- 
alent widths, W r = 0.65T0.10 and 0.58±0.10A respectively 
at z g . Using the empirical rel ationship between E(B- V) and 
W(Na i) reported recently bv lPoznanski et"al] (|2012h we es- 
timate E(B-V) = 0.31±o:i2 and 0A7±%]£ from the measured 
equivalent widths. The E(B-V) we measure from the SED 
fitting is higher than these values but consistent within un- 
certainties. 

As the doublet ratio is close to 1 and the resolution of 
our spectrum is low it is most likely that the lines are satu- 
rated. A very conservative lower limit of log A(Na i)^ 12.78 
is obtained assuming the optically thin case. Consequently 
we derive A r (H i)^ 10 21 cm~ 2 from the known correla- 
tion between A(Na i) and A(H i) found in our Galaxy 
(|Ferlet et alj|l985l : IWakker fc Mathisll2000r ). 



2.4 Diffuse Interstellar Bands 

All the optical observations presented above suggest that the 
gas along the line of sight has A(H i) 10 21 cm~ 2 , with near 
solar metallicity and high dust content. In such cases one can 
expect to detect absorption fr om differen t molecular species 
and diffuse interstellar bands (|Sarrell2006h . In Fig. [2] we have 
marked the expected positions of DIBs using vertical dashed 
lines. A clear absorption is detected at the expected position 
of the A r = 6284 A DIB. Using a Gaussian fit we measure a 
rest equivalent width of 1.46T0.21 A and FWHM = 14±2 A. 
The measured FWHM is 1.4 times larger than what is seen 
in nearby starburst galaxies (|Heckman fc Lehnertll200ol ) and 
roughly a factor 2.2 high er than what is seen in the Galactic 
ISM towards HD 204827 (|Hobbs et alj|2008l ). No other DIB 
absorption is detected at more than the 3cr significance level. 
We place a 3er upper limit of 0.66 A for the rest equivalent 
width of the A5780 DIB feature. 

In the local universe DIB obser vations are avail- 
able for the ISM of M agellani c clouds dWelty et all 120061 ; 
ICox et al.1 12006| I2007T ). M33 (ICordiner et al.1 |200ST) M31 
( Cordiner etldl I2011T ), NGC 1448 dSollerman et al.1 12005T) 
and some nearby starburst galaxies (|Heckman fc Lehnertl 
l2000h . In Fig.E]we plot W r (A6284) against E(B-V) found in 
these cases. The strength of A6284 is affected by the strength 
of the radiation field and metallicity. Low equivalent width 
is explained by stronger radiation field and/or lower metal- 
licity It is interesting to note the metallicity measured in 
the present case is close to that measured in M31 an d star- 
burst galaxies studied bv lHeckman fc Lehnertl (|2000h . How- 
ever, the equivalent width of the A6284 DIB is found to 
be higher than t hat p redicted by the correlation found by 
IFriedman" et al.1 <|201ll ). This may imply slightly low ambi- 
ent radiation field along the QSO sight line. Alternatively, 
if we assume that the physical conditions in the galaxy 
studied here are similar to that of the Galactic ISM then 
we get log JV(H i) = 21.76T0.25 and E(B-V) = 1.30T0.20 
fro m the observed Vrv(A6 284) and the correlations found 
by IFriedman et al.l (|201ll ). Interestingly the N(R i) value 
inferred here matches well, albeit with large uncertainties, 
with that w e infer from the Av— A(H i) relationship of 
iBohlin et all (|l97Sf ). In addition, the inferred E(B-V) may 
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and not 3.1 as we have assumed before. 
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Figure 3. Comparison of E(B-V) and 1VV(A6284) in differ- 
ent extra-galactic environments. The solid and dashed lines are 
the best fi tted relationship and its 1 a error found in the 
milky way l|Friedman et al,1l201lh . Filled tria ngles, open square s 
and open circles are for Magellanic clouds jWeltv et al J [2006ft . 
nearby starburst gala xies ( Hc ckman fc Lehnertl l2000h and M31 
iCordiner et al.ll201lft respectively. The data for the present sys- 
tem is marked with a star. 



This is similar to what h as been inferred fo r the host galax- 
ies of high-z supernovae jWanget al.ll2006T l. Therefore, it is 
important to have independent measurements of iV(H i) and 
Rv for the present case. 

In all the local measurements discussed above 
W r (A5780) is found to be smaller than jj^ 1X628 4). If we use 
the correlations found bv lFriedman et alj (|201lT ), we obtain 
VF r (A5780) = 0.61 ±0.08 A. We show the expected profile in 
the bottom panel of Fig. [2] From the residual plot shown in 
the top panel of Fig.[5]we can say that the observed spectrum 
is consistent with A5780 DIB featu re having th e pred icted 
equivalent width. From Table 3 of IWeltv etail (|200ot ). we 
find the average value of the W r (A6284)/W r (A5780) ratio to 
be 2.9T0.2 in the case of the Magellanic clouds. Based on 
this we obtain VK(A5780) = 0.5 ±0.1 A. Detecting this fea- 
ture and other DIB features in a high signal-to-noise ratio 
spectrum will allow us to probe the physical state of the gas 
in more detail. 

This is only the third detectio n of DIBs due to interven- 
ing absorbers in Q SO spectra (see I Junkkarinen et alj|2004l ; 
lEUison et alJlioOa ). In previous cases the A5780 DIB was 
detected without clear detection of the A6284 DIB. This 
is very much contrary to what has been seen in the local 
univer se (apart from the rare sight line towards Sk 143, as 
seen m IWeltv et al.ll2006t ) and in the present case. lYork et al.l 
(|2006ft . by analogy with Skl43, attributed these unusual line 
ratios to the ISM being more protected from the ambient UV 
radiation field. 

The present NTT spectrum does not allow us to search 
for molecular absorption. However, based on known corre- 
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Figure 4. Left panel: Velocity plot showing the Ha emission ex- 
tracted at different locations in the galaxy. Top right: 2D spectra 
in the region of Ho and [N if] emission lines after subtracting 
the QSO light. The horizontal lines show the location of the QSO 
trace. Bottom right: the rotation curve obtained from the Ho line. 
The star is the velocity of the H I gas seen in 21-cm absorption 
towards the QSO (see Section |3,1| I, The zero of the velocity scale 
is set at 2 a b s = 0.0791 and zero of the spatial scale is set at the 
position of the QSO. 



lations in the Galactic ISM we expect log iV(CH)^13.8 , 
log iV(CH + )5;13.5 and log A(CN)^12.4 jWeltv et alj|2006h . 
Therefore high resolution follow-up spectroscopy of this 
source could yield very good insights into the ISM of this 
external galaxy. 



2.5 Ha rotation curve 

In this sub-section we study the large scale kinematics of 
the emission line gas using the NTT long-slit spectrum. In 
the left panel of Fig U we show the Ha emission observed in 
different sub-apertures together with the best fitted Gaus- 
sian profiles. The measured FWHM is close to our spectral 
resolution suggesting low velocity dispersion along the line 
of sight. The 2D spectrum around the Ha range after sub- 
tracting the QSO trace is shown in the right panel. While 
the rotation is apparent we also notice that the blue side is 
brighter suggesting additional star formation activity. Coin- 
cidentally the QSO sightline is very close to this star form- 
ing region. The angular separation between the center of 
the QSO trace and the center of the nearest star forming re- 
gion (based on the peak of the Ha emission) is 0.48 arc sec. 
This corresponds to a physical separation of 0.7 kpc. The 
different values of Ay measured toward the Ha emitting re- 
gion and the QSO reflect roughly a factor of 2 change in 
Ay within a projected separation of 0.7 kpc. Therefore, the 
projected distance (and associated dust extinction) is large 
enough so that the gas along the QSO sightline may not be 
influenced by the star-forming region. The QSO sight line is 
at an impact parameter of 4 kpc from the galactic center. 
The rotation curve plotted in Fig. [4] suggests an asymptotic 
circular velocity of ~ 125 km s _1 . This circular velocity for 
a typical radius of 10 kpc corresponds to a dynamical mass 
of ~ 4 x 10 10 M© after correction for the inclination angle 
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Figure 5. Optical depth profile of the 21-cm absorption towards 
J1639+1127 detected in the GMRT spectrum. Top panel shows 
the cumulative optical depth profile. The open squares (and ver- 
tical dashed lines in the bottom panel) mark the velocities at 
which the integrated optical depth is 5%, 50% and 95% of the 
total value. The mean optical depth and associated RMS in the 
absorption free regions are shown with solid and dotted horizon- 
tal lines respectively. The velocity scale is defined with respect to 
z = 0.079098. 

of 63 deg. This suggests that the host galaxy is a low mass 
disk galaxy. 



3 21-CM ABSORPTION AND PARSEC SCALE 
STRUCTURES 

3.1 GMRT 21-cm spectrum 

As the background QSO is radio-loud (with 1.4 MHz flux 
density of ~ 164 mjy) we get an opportunity to study the 
property of the cold H I gas along this dusty sightline. The 
source is unresolved in the FIRST (Faint Images of Radio 
Sky at Twenty-centimeter) image with a deconvolved size of 
1.32" xO. 82". As expected, a very strong 21-cm absorption 
was detected by Giant Metrewave Radio Telescope (GMRT) 
observations on 02 July 2010 (2 MHz bandwidth split into 
128 channels having a velocity resolution of ~3.5 km s _1 per 
channel). We obtained a higher resolution (2.08 MHz band- 
width split into 512 channels having a velocity resolution of 
0.93 km s _1 per channel) GMRT spectrum on 14 and 16 July 
2011. All these data were reduced using the National Radio 
Astronomy Observatory (NRAO) Astronomical Image Pro- 
cessing System (AIPS) following the standard procedures. 
The radio source is unresolved in our GMRT image (with 
a spatial resolution of 2.5" x 2.1") with a peak flux density 
of 155 mjy/beam. The 21-cm optical depth profile from the 
high resolution spectrum is shown in Fig. [5] The integrated 
21-cm optical depth is found to be, T =15.7±0.13 km s _1 , 
with 90% of it is within a velocity range of 28 km s~ (see 
top panel in Fig. [5}. The observed T translates to 7V(H i) 
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Figure 6. VLBA image showing J1639+1127 resolved into 3 
Gaussian components. Crosses mark the locations of components 
A, B and C from top to bottom. 21-cm H I spectra at the peak 
of two strongest components, A and C, separated by 89 pc, are 
also shown. Absorption profiles are clearly different. The RMS 
in the image is 8.425 xl0 — 5 Jy/beam and the contour levels are 
4.0xl0~ 4 x (-2,-1,1,2,4,8,16,32, ...) Jy/beam. The restoring beam 
(0.01271" xO. 00533" with a position angle is -3.18 deg) is shown 
as an ellipse at the bottom of the figure. Image center is at RA 
= 16 39 56.3624 and DEC = +11 27 58.6568. 

= 2.9 x 10 19 (T S // C ) cm" 2 with f c being the covering fac- 
tor, the fraction of the background radio emission covered 
by the absorbing gas. The probable values of A(H i) dis- 
cussed above based on known correlations suggest that the 
harmonic mean spin temperature of the H I gas is less than 
or equal to 100 K as seen in the Galactic ISM. 

Another interesting observation is that the H I gas that 
produces the 21-cm absorption is redshifted with respect to 
the peak rotational velocity of the Ha emitting gas at the 
same location by up to 25 km s _1 (see Fig. [4]). This once 
again confirms that the absorbing gas and Ha emitting gas 
are not co-spatial. 



3.2 Spatially resolved VLBA spectroscopy 

The VLBA observations of J1639+1127 were carried out on 
July 27 and October 1, 2011. The total and on-source ob- 
serving times were 10 and 6.3 hrs respectively. One 8 MHz 
baseband channel pair was used in these observations, with 
right- and left-hand circular polarization sampled at 2 bits, 
and centered at the frequency of the redshifted H I 21-cm 
line. The data were correlated with 4096 channels (chan- 
nel width of 0.45 km s _1 ) and 2 sec correlator integration 
time. We followed the standard data reduc tion procedures 
for reducing the VLBA dat a (see for example iMomiian et al.l 
|2002| ; ISrianand et"al]|2012h . The radio continuum emission 
is resolved into multiple components at mas scales (Fig. |BJ). 
We use the multiple Gaussian fits to the VLBA images and 
found 3 distinct components. The results are summarised in 
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Figure 7. Bottom three panels show the 21-cm optical depth 
profiles observed in, respectively, the VLBA spectra of compo- 
nents A, C, and the GMRT spectrum. The second panel from the 
top shows the optical depth difference observed between the two 
VLBA sightlines. The top panel shows the fractional change in 
the optical depth between the two sight lines. The velocity scale 
is defined with respect to z = 0.079098. Vertical dashed lines are 
as defined in Fig. [5] and the dotted lines mark the velocities of 
narrow components seen towards C. 



Table 2. Results of VLBA observations 



Source Flux density Deconvolved 

Peak Total angular size PA J T(v)dv 

mjy/beam mjy (mas 2 ) (deg) kms — 1 

A 26.8 57.4 9.9x6.7 +31.6 17.6+1.6 

B 1.7 2.1 8.8x0.0t +107.2 

C 36.8 56.8 5.9x4.1 +62.0 37.6±2.9 



tion components superposed on a smooth absorption profile. 
This confirms the patchy distribution of the absorbing gas 
at parsec scales. 

To probe this further, we compare in Fig. [7] the optical 
depth profiles seen in the GMRT and the VLBA observa- 
tions. Fig. [JJ also shows the difference in optical depth (i.e 
At) between the lines of sight towards A and C as a func- 
tion of velocity. While broad smooth absorption profiles with 
nearly identical optical depths are seen, narrow absorption 
components (identified with vertical dotted lines in Fig. [TJl 
are clearly seen only towards C. The top panel of Fig. [7] 
shows the fractional difference in r between the two sight 
lines compared to tc ■ It is clear that on average tq ~ 5 x ta 
in the velocity range -5 to 10 km s~ , where the optical depth 
towards C is dominated by several narrow components. This 
basically confirms the large variations in H I optical depth 
over a length scale of 89 pc for narrow H I components. 

3.3 Tiny clouds and parsec scale structure 

As the deconvolved sizes of the VLBA components corre- 
spond to a projected size of ~ 8 pc the next question we 
wish to address is whether there is any optical depth vari- 
ation over this length scale. We do this by comparing the 
optical depth profiles seen towards A and C with the total 
optical depth profile seen in the GMRT spectrum. In general 
we can write, 

TGMRT(v) = fATA.(v) + fcTc(v) + Td(v). (2) 

Here, fA and fc are the covering factors of the gas com- 
ponent having optical depth ta(v) and tc(v) respectively. 
td(v) is the optical depth towards components not seen in 
the VLBA image, fc will be 0.37 if the absorbing gas cov- 
ers all the radio emission seen in the VLBA image from 
the component C (i.e fc is ratio of flux density of compo- 
nent C to the total flux density measured in the GMRT 
image). Therefore, if we assume a plane parallel gas slab cov- 
ering only the emission from component C then we expect 
tgmrt(v) = 0.37tc(v). If we assume ta(v) = tc(v), then we 
expect tgmrt(v) = 0.74rc(v). If Ta(v) tc(v) then we ex- 
pect 0.37 ^ (tgmrt(v)/tc(v)) < 0.74. On the contrary if 
Ta(v) > tc(v) then tgmrt(v)/tc(v) ^ 0.74. Note this condi- 
tion is also obtained when additional components are present 
towards the broad diffuse emission resolved out in the VLBA 
image (i.e td 0). This discussion clearly suggests that by 
looking at the ratio of the optical depth observed by GMRT 
and VLBA in one of the components we will be able to 
draw broad conclusions about the optical depth variability 
at small scales. 

First we focus on the most interesting narrow compo- 
nent at v = km s _1 . This component is distinctly visible in 
C and evident even in the GMRT spectrum however clearly 
absent towards A (as suggested by At ~ tc)- The peak op- 
tical depth of this component is 2.5T0.1 towards C. Based 
on the flux densities in the VLBA image, if the absorbing gas 
covers only C then we expect the peak optical depth at v ~ 
km s _1 in the GMRT spectrum to be ~ 0.93±0.05. However, 
the measured value in the GMRT spectrum is 0.70 ± 0.01, 
suggesting that, if the absorption component is a plane par- 
allel slab with uniform JV(H i), then it would cover only 75% 
of the continuum emission from component C. The observed 
parameters, on the contrary simply imply the presence of a 



t An angular extent of indicates a source size much smaller 
than the synthesized beam. 

Table [2] Only 78% of the flux density measured with the 
GMRT is recovered with the VLBA. Components A and 
C, that are separated by 55.8 mas (implying a linear line 
of sight separation of 89 pc at the redshift of the galaxy), 
have nearly equal flux densities and contribute to 98% of 
the total flux density seen in the VLBA image. The peak 
flux densities of these components are 47% and 65% of their 
respective total flux densities, suggesting structures at few 
mas scales in the radio emission. 

The 21-cm absorption toward the peak emission of A 
and C extracted from our VLBA observations (smoothed to 
a channel width of 0.9 km s~ to match the GMRT spectrum) 
are shown in Fig [5] The integrated 21-cm optical depths, 
given in the last column of the Tabled are more than that 
measured in our GMRT spectrum and differ by a factor of 
2 within the 89 pc probed by the two sightlines. Unlike the 
GMRT spectrum, that is very smooth, the 21-cm spectrum 
towards component C shows the presence of narrow absorp- 
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Figure 8. The ratio of tqmrt and tc as a function of relative 
velocity. The velocity scale and vertical lines are as in Fig. [7] 



strong optical depth gradient within the deconvolved size of 
component C (i.e ^6 mas). This angular scale corresponds 
to a transverse size < 8.8 pc. As this component is distinctly 
visible, we fit multiple Gaussian to tc- While Gaussian fit 
to other components may be unphysical, simultaneous fits 
allow us to get a realistic estimate of FWHM (1.75 km s _1 ) 
and the peak optical depth (2.5T0.1) for the narrow compo- 
nent. We get an upper limit of 66 K for the kinetic tempera- 
ture of the gas from the FWHM. Using this as an indicator of 
spin temperature (Tg), we get iV(H i) = 5.4 x 10 20 (T S /66K) 
cm~ 2 . If we approximate the absorbing cloud as a sphere 
then we get a particle density of ~40 cm~ 3 using iV(H i) 
and the transverse size discussed above. Physical conditions 
in this component are typical of diffuse molecular clouds. 

To explore this further we plot in Fig. [8] the ratio of 
tgmrt and rc in the velocity range that contains 90% of the 
integrated optical depth. The expected ratio when a parallel 
slab of gas with constant r covers only the component C, all 
the VLBA components and all the flux seen by GMRT are 
shown by horizontal lines at 0.37, 0.74 and 1 respectively. We 
use this plot to discuss the extent of the gas that produces 
absorption in different velocity ranges. 

First we consider the velocity range 3 ^ w(km s" 1 )^ 9 
where absorption is seen towards both A and C with tc ~ 
5 x ta (see top panel in Fig. [7|). If the gas towards C and 
A covers all the radio emission from these components then 
we expect tgmht/tc to be 0.44 (i.e. 0.37 x (1.0 + 0.2)). 
From Fig. [8] we see the observed ratio is consistent with 
0.38 in this range. We notice that the optical depth error 
in each channel in the VLBA spectrum towards C is ~0.1. 
For the mean optical depth measured in this velocity range 
this translate to an error of 0.03 in the ratio. Therefore, in 
each channel we find the ratio to be lower than 0.44 by 2a 
level. As this happens over 6 channels we see the difference 
to be significant at the 4.8 a level. This difference can be 
understood if some of the narrow components seen towards 
C have projected sizes less than 8.8 pc (or the presence of 
strong opacity gradient within this scale) as seen in the case 
of the narrow component at v ~ km s~ x . 

In the velocity ranges 10 ^ w(km s -1 )^ 16 and — 10 ^ 
w(km s -1 )^ —4 from Fig[7] we find At ~ 0. Therefore we 
expect the ratio of tgmrt and tc to be 0.74. In the case 



of the first velocity range the ratio is above 0.74. In the 
framework presented above, this could either mean ta 
tc or contribution to tgmrt from the diffuse component 
resolved in the VLBA image. If we use the average value 
(1.2) of the measured ratio in this velocity range we get 
ta ~ 2.1 x tc ■ In the second velocity interval the ratio is 
found to be between 0.36 and 0.74. This is consistent with 
t a ^ tc- Using the mean value of the ratio (0.65) in this 
velocity range we find r c ~ \Ata- These are very different 
from tc = ta we see towards the peak emission in A and 
C. Therefore, we can conclude that there are strong opacity 
gradients at the spatial scale of ~ 10 pc even in the gas that 
produces absorption in this velocity range. 

Finally in the velocity range -18 to -10 km s~ , the 
absorption profile is smooth both in the VLBA spectrum 
towards C and the GMRT spectrum. Therefore we smooth 
the spectra to a 4 km s" 1 resolution to increase the signal to 
noise ratio. We notice that tgmrt/tc in this velocity range 
is 0.67±0.03. This just differs from the expected value of 
0.74 by 2.3a level. Within the observational uncertainty the 
optical depth seen towards A and C are nearly equal as sug- 
gested by the difference in r plotted in Fig. [7] Increasing 
the signal to noise of our VLBA spectra will help detecting 
minor optical depth differences even in this smooth compo- 
nent. 

In summary, we find the 21-cm absorption has both nar- 
row and broad components. The narrow components show 
optical depth variations by a factor 1.4 to 10 over a 89 pc 
scale. In one of the narrow components we show that optical 
depth variations are present over a 8 pc scale. On the con- 
trary the absorption seen in the broad wings are consistent 
with nearly similar optical depths (i.e at 2.5<r level) within 
measurement uncertainty. Thus a simple picture of the sys- 
tem could be that of several cold, dense and small clouds 
with characteristic size ^ 10 pc embedded in a smooth dif- 
fuse H I medium covering several tens of parsec. This sce- 
nario is_similax_to_the "bjobby sheets" seen in the Galactic 
ISM (|Heiles fc Trolandll2003l ). 



4 DISCUSSION 

We report the detection of diffuse interstellar bands 
and 21-cm absorption from the z a b s = 0.079 galaxy 
SDSS J163956.38+112802.1 towards the z om = 0.993 QSO 
J163956. 35+112758.7. The QSO line of sight in this case is 
passing through the disk of a relatively low mass (~4 x 10 10 
Mo), near solar metallicity spiral galaxy at an impact pa- 
rameter of ~4 kpc. The region with highest star formation 
surface density, inferred from Ha emission, in this galaxy is 
located within a projected distance of 0.7 kpc from the QSO 
sightline. 

The QSO is highly reddened with E(B-V)~0.83±0.11 
(estimated from the SED fitting) along the line of sight. 
This suggests that the line of sight is passing through a 
translucent interstellar medium. The measured E(B-V), rest 
equivalent widths of Na D lines and surface density of star 
formation favor N(H i) ^ 10 21 cm -2 along this sight line. 

We report the detection of a A6284 DIB feature in the 
QSO spectra at the redshift of the galaxy. The strength of 
this feature per unit reddening is found to be higher than 
what is seen in the Milkyway. This could be related to dif- 
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ferences in the DIB carrier and/or to the ambient radiation 
field. We also find the FWHM of this feature to be broader 
than what is typically seen in our Galaxy but close to what 
is seen in starburst galaxies. Our spectrum is consistent with 
the presence of the A5780 DIB with a rest equivalent width 
predicted from the ob served A6284 DIB featu re and the cor- 
relation reported by iFriedman et all (|201ll ). Due to poor 
SNR, no other DIB feature is detected in our spectrum. In 
order to probe the physical conditions and the nature of the 
DIB carriers it would be important to obtain high signal-to- 
noise ratio spectrum of the QSO. 

Strong 21-cm absorption from the galaxy, spread over 28 
km s _1 , is detected in our GMRT spectrum. Our follow-up 
VLBA observations reveal the presence of two strong radio 
emitting components, separated by 89 pc, that are used to 
understand the spatial variations of the 21-cm optical depth. 
First we notice that the integrated optical depths towards 
these two components differ by a factor of 2. We show that 
this difference is mainly dominated by narrow components 
that are seen only towards one of the VLBA component (i.e 
component C). In these narrow components the H I opacity 
seems to change over ~ 10 pc scales. For one of these narrow 
components we show the kinetic temperature is ^ 66 K. As 
we expect dust to be mainly associated with low tempera- 
ture gas, if the other narrow components are also cold, then 
we can speculate that the optical emission from the QSO is 
associated with the VLBA component C. 

In high-z 21-cm absorption line studies of DLAs and 
Mg II systems, and in the absence of VLBI spectroscopy, one 
uses the 'core' fraction measured from the mas scale contin- 
uum images to co rrect the optical depths for partial coverage 
(see for example, Kanekar et al]|2009l ; ISrianand et al-lfeoi^ 
iGupta et al.112012 ). In case of DLAs, this allows one to mea- 
sure the spin temperature of the H I gas. There are three 
issues that affect this approach: (i) Unambiguous identifica- 
tion of the 'core'. If the milliarcsecond morphology of the 
radio emission is not simple, then VLBI observations are 
needed to identify the flat-spectrum 'core' component, (ii) 
The fact that the radio sight line towards the 'core' may 
still be tracing a larger gas volume with respect to the op- 
tical sight line, (iii) The implicit assumption that a single 
absorbing cloud covers the 'core' i.e. a single covering factor 

(/c). 

The third issue is mainly related to the patchiness of the 
absorbing gas on scales of tens of pc, which is typically the 
resolution provided by VLBI observations, and can only be 
addressed via VLBI spectroscopy. In the case of J1639+1127 
discussed here it is plausible to associate component C with 
the quasar sight line because it shows clear signatures of the 
presence of a large amount of cold gas. Under the assump- 
tion of a single cloud covering C we would have estimated f c 
as 0.37 (57/155). However, if some of the absorption com- 
ponents were to cover only the peak of C, then f c would be 
off by as much as a factor of 1.5 (56.8/36.8) from the above 
value. Furthermore, the covering factor estimate could be 
dramatically off if the core is associated with component B, 
resulting in f c of 0.01, or with a component not detected in 
our VLBA image. This is a possibility considering that the 
radio continuum morphology of J1639+1127 at high angular 
resolution (se e Fig[6l) resemb les that of a compact symmetric 
object (CSO; IConwavll2002l ). Here, the two dominant radio 
sources, A and C, would be the two radio lobes of the CSO. 



Therefore, mult i- frequency VLBI observations are needed to 
detect and/or identify the core component to address this 
issue. 

The case of J1639+1127 also demonstrates that asso- 
ciating the entire absorption detected at arcsecond scales 
through GMRT observations with either A or C component 
would lead to f c estimates that are at least off by a fac- 
tor of 2. In general, this means that even when the core 
component is accurately identified, Ts measurements from 
DLA studies can have large errors especially if the back- 
ground quasar has significant radio structure at the scales 
where the absorbing gas has strong optical depth gradients. 
Therefore one requires more high-z measurements towards 
strong core dominated sources to address the issues related 
to the measurements of Ts- 

Based on the measured E(B-V), the equivalent widths 
of A6284 DIB and Na D lines and various correlations found 
in our Galaxy we predict appreciable column densities of 
molecular species like CH, CH + and CN along this line of 
sight. Detecting these species will allow us to probe the phys- 
ical and chemical state of this translucent gas in detail. In 
addition, this is a good target to search for OH lines and 
molecular lines in radio/mm wavebands. 
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